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We study through numerical simulation the spherical collapse of isothermal gas in Newtonian 
gravity. We observe a critical behavior which occurs at the threshold of gravitational instability 
leading to core formation. This was predicted in a previous work by two of the present authors. We 
describe it in detail in this work. For a given initial density profile, we find a critical temperature T* , 
which is of the same order as the virial temperature of the initial configuration. For the exact critical 
temperature, the collapse converges to a self-similar form, the first member in Hunter's family of self- 
similar solutions. For a temperature close to T* , the collapse first approaches this critical solution. 
Later on, in the supercritical case (T < T*), the collapse converges to another self-similar solution, 
which is called the Larson- Penston solution. In the subcritical case (T > T*), the gas bounces and 

fC) , disperses to infinity. We find two scaling laws with respect to [T — T*|: one for the collapsed mass 

in the supercritical case and the other, which was not predicted before, for the maximum density 
reached before dispersal in the subcritical case. The value of the critical exponent is measured to 
be ~ 0.11 in the supercritical case, which agrees well with the predicted value ~ 0.10567. These 
j.^ ■ critical properties are quite similar to those observed in the collapse of a radiation fluid in general 

' relativity. We study the response of the system to temperature fluctuation and discuss astrophysical 

^5 | implications for the interstellar medium structure and for the star formation process. Newtonian 

critical behavior is important not only because it provides a simple model for general relativity but 
also because it is relevant for astrophysical systems such as molecular clouds. 

<n : 

PACS numbers: 04.40.-b, 64.60.Ak, 97.10.Xq 

m ■ 
> ■ 

t^" ! I. INTRODUCTION 

O ■ 

The Newtonian isothermal gas is a relevant description for important astrophysical systems. It is indeed a good 
7— I ' description for cold molecular clouds found in the galactic interstellar medium. In these clouds the typical cooling time 
is usually much shorter than the dynamical time; as a result they are in thermal equilibrium with whatever steady 
heating source is present, for example the blackbody radiation of the cosmic microwave background, extragalactic 
ultraviolet radiation, or local stars. The molecular clouds can be described as isothermal. Since they are the prime 
CP location for star formation, their physics determines the whole star formation history in the galaxy. A detailed 
5^ | understanding of the physics of these systems will lead to predictions for the properties of star formation. 
&)[)■ There are two possible theoretical approaches for the self-gravitating Newtonian gas: statistical physics and dy- 
namics. The traditional statistical physics approach deals with equilibrium configurations of the system and gave the 
•rH , first description of Jeans gravitational instability which is a fundamental behavior in the physics of self-gravitating 
systems (Chandrasekhar [|, Padmanabhan 0). Using the modern tools of statistical field theory and renormalization 
groups, de Vega, Sanchez and Combes Q and Semelin et al. give a theoretical basis to the observed self-similar 
properties of the gas. These statistical methods are useful, generally speaking, to describe systems undergoing second 
order phase transition and exhibiting critical behavior. Investigating the critical behavior of the gas is also the subject 
of the present work, where we choose to describe the system through dynamics rather than statistical physics and 
through the renormalization group. 

The basic behavior of a clump of isothermal gas is the following. In a sufficiently dense and/or cold configuration, 
the gas will collapse due to its self-gravity, while in a diffuse and/or hot one it will expand due to its internal pressure 
forces. A critical behavior appears at the transition between these two regimes. 

In fact, Choptuik Q showed the existence of a so-called critical phenomenon in the spherical collapse of a massless 
scalar field by general rclativistic numerical simulations. The simulations were performed for a one-parameter family 
of initial data sets. He found a self-similar solution with critical properties at the threshold of black hole formation. 
In this first case, the self-similarity holds for a discrete set of transformations only (periodicity in logarithmic time). 
He showed that the mass of the formed black hole follows a scaling law, and that the phenomenon is universal in 
the sense that it does not depend on the one-parameter family of initial data sets. Soon after that, similar critical 
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phenomena were observed in the spherical collapse of radiation fluid by Evans and Coleman [6j. In this second 
case, self-similarity holds for a continuous set of transformations. Koike, Hara and Adachi gave a clear picture 
of these phenomena by using the renormalization grou p ap proach. Recently, critical phenomena have been found in 
perfect fluids with more general equation of state |8L I9L llOj] and in a variety of systems. Astrophysical implications of 
the critical phenomena have been studied in the scenario of primordial black hole formation by Yokoyama and 
Niemeyer and Jedamzik |l2j . See [Tsl ] for a recent review of critical phenomena in gravitational collapse. 

The solution with critical properties is only one of several known self-similar solutions. See Carr and Coley [T^| 
for a recent review of self-similar solutions in general relativity. It was conjectured by Carr |15| that a generic 
solution gradually evolves toward a self-similar form, which suggest that one of the self-similar solutions may be an 
attractor solution. This is called the self-similarity hypothesis. In Newtonian gravity, a self-similar solution, the 
so-called Larson-Penston solution 0, is known to be an attractor solution in the spherical collapse of isothermal 
gas 0, 0, |2jj. In general relativity, Harada and Maeda [2l] showed by numerical simulations that the general 
relativistic counterpart of the Larson-Penston solution, which was discovered by Ori and Piran j2^ , is an attractor 
solution in the spherical collapse of a perfect fluid with P = ke (0 < k < 0.03), where P and e are the pressure 
and energy density, respectively. This convergence is quite important in general relativity because it provides a 
strong counterexample against the cosmic censorship conjecture in spherical collapse 0, |22, |U 0] • See [25| for a 
recent review of cosmic censorship and related topics. In the context of the self-similarity hypothesis, the stability 
criterion for general relativistic self-similar solutions against kink mode was recently obtained by Harada |26| , through 
a generalization of the Newtonian analysis j2]j • 

The present paper shows that critical phenomena exist in Newtonian gravity also. We concentrate on a spherically 
symmetric system of isothermal gas in Newtonian gravity. The existence of a critical phenomenon was predicted by 
Harada and Maeda [2l| and Maeda and Harada [23 • In these papers, based on the renormalization group analysis, 
it was predicted that a critical phenomenon is associated with the first member of Hunter's family of continuous 
self-similar solutions |2^. The criticality appears in the scaling law for the collapsed mass, which is satisfied for 
continuous values of the mass. The value of the critical exponent was predicted to be ~ 0.10567. The critical 
phenomenon observed in this paper agrees very well with the above predictions, and gives a precise meaning to the 
collapsed mass in Newtonian gravity. We have also observed that the scaling law holds over many orders of magnitude 
in terms of the parameter. 

The organization of this paper is as follows. In Sec. II, we apply the renormalization group analysis to the 
Newtonian isothermal gas system and draw the expected properties for the critical phenomenon. In Sec. Ill, we 
present the results from numerical simulations of Newtonian collapse. In Sec. IV we study the response of the system 
to temperature fluctuations and discuss astrophysical applications. In Sec. V, we summarize the paper. 



II. RENORMALIZATION GROUP ANALYSIS 



In this section, we review the renormalization group analysis of the critical behavior in a spherical system of 
Newtonian isothermal gas. Part of the following discussion is based on Maeda and Harada psj . 

The dynamics of a spherically symmetric isothermal gas system is described in the following set of equations: 

9 , „ 19.99. 9 dp Gm 
dm dm 

~dt +V ^ = ^ (2 - 3) 
dm 9 

— = 4nr 2 p, (2.4) 

where p = p(t,r), v — v(t,r), m — m(t,r), c s and G denote the density, radial velocity, total mass inside the radial 
coordinate, sound speed, and gravitational constant, respectively. The sound speed c s is related to the temperature 
T of the gas as 

c s = (fcT/m) 1 / 2 , . (2.5) 

where k and m are the Boltzmann constant and mass of each particle, respectively. 

In order to apply the renormalization group approach, it is convenient to redefine the physical quantities as follows: 

r = — , v = — , m = — . (2-6) 

c s c s c„ 
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In terms of these new quantities, the basic equations reduce to the following form: 

^W + ^trVl^fp^ , (2.8) 
flF + *»r = ' < 2 - 9 > 

— = 47rf 2 p- 2.10) 
or 

We can see that the above hatted system does not contain any dependence on the temperature (or sound speed), 
which simplifies the analysis. 

The initial conditions for the bare system are specified by the density and velocity profiles at the initial time 
t = t\ < (t — is the time when the density singularity occurs in the center or, more physically, the time of core 
formation). Suppose the initial density and velocity profiles are given by functions p\{r) and v- 1 {r) 1 

p{t h r) = Pi (r), v(t u r)=Vi(r), (2.11) 

respectively. Using the relations of the hatted quantities to the bare quantities, we specify the initial conditions for 
the hatted system by 

p(t u r) = p^f), fi(ii,f) = ^M. (2.12) 

Cs 

Therefore, the parameter-dependent system (I2.1|) - (|2.4|) with unique initial conditions is equivalent to the parameter- 
independent system (|2.7|I - (|2.1U|) with a one-parameter family of initial conditions. The parameter is c s or T. 

Although we avoid to fully repeat the analysis by Koike, Hara and Adachi Q, we need to review it briefly. We 
introduce the variables 

T = -ln(-t), x = ln—. (2.13) 

r is the scaling variable in the renormalization group transformation. Suppose a self-similar solution H ss (x) exists, 
with a unique relevant (unstable) mode, which turns out to be the critical solution. This self-similar solution can be 
considered as a fixed point of the renormalization group transformation. The uniqueness of the unstable mode implies 
that the fixed point H ss has a stable manifold of codimension one. This manifold is referred to as a critical surface. Any 
one-parameter family {Hr p \(x)\p € R} of initial data sets with parameter p generically has an intersection H* with 
the critical surface of the fixed point H ss . H* will be driven to -ff ss under the renormalization group transformation, 
and hence H* is the initial data set with the critical value p* . We consider an initial data Hi{x) in the one-parameter 
family, close to H* , i.e., 

h(0,x)=Hi(x)=H*(x) + eF(x), (2.14) 

where h(r, x) is a solution of the partially differential equations and e = p — p* . Then, for large r , we have 

h( T0 ,x) ~ H ss (x) + ee KT °F iel (x), (2.15) 

where K and F Te \(x) are the eigenvalue and eigenfunction of the relevant mode, respectively. Now we choose tq so 
that the first term and the second term are of the same order, i.e., 

ee KT0 = 0(1). (2.16) 

At r , the solution is deviating from the critical solution. If the collapse bounces to expansion, which will be the case 
for a subcritical collapse, the above condition obviously holds for the velocity field at the moment of the bounce. On 
the other hand, if the finite mass near the center shrinks considerably faster than for the critical collapse, which will 
be the case for the supercritical case, we have a tq for which the above condition also holds. Since the deviation from 
the critical collapse grows at x = 0(1), the relation f = e x ~ T implies that the radius of the core collapse or bounce is 
given by f = O(e~ T0 ) = 0(e 7 ), where 7 = l/(Re«). 

In order to see clearly that self-similar solutions exist in the present hatted system, we introduce the following 
quantities: 

U = v, P = Attt 2 P , M=— ^. (2.17) 
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Then, the equations transform to 

-P + (1 + zU)P' + zPU' = 0, (2.18) 

(UP + UP) -{U + zU 2 + z)P' - (1 + 2zU)PU' - 2zP + MPz 2 = 0, (2.19) 

M - M' - M + PU = 0, (2.20) 

-zM' = P, (2.21) 

where the dot and prime denote the partial derivatives with respect to r and x, respectively, and z = e~~ x . The 
self-similar solutions are obtained by numerical integration of the ordinary differential equations which follow from 
the above equations with the assumption P = P(x), U = U(x) and M = M(x). If we assume the existence of the 
critical solution, using the relation 

c 3 

m = ^M(x)e~ x r, (2.22) 

we obtain the scaling law for the collapsed mass or the "core" mass for the supercritical collapse, 

m corc oc \p- p*p. (2.23) 



For the subcritical case, using the relation 



P(aO 



r~\ (2.24) 



47T 

we obtain the scaling law for the maximum density, 

p max a|p-pT 2 T. (2.25) 

In fact, for Eqs. (|2.18() - (|2.21() . there is a discrete set of self-similar solutions with analyticity at the center and at 
the sonic point, such as the homogeneous free-fall, Larson-Penston, Hunter (a), (b), (c) and (d) solutions, and so on. 
Larson-Penston solution describes the coherent collapse of the cloud for t < 0. This solution can be extended beyond 
t = to a late-time solution in which a finite mass collapses to the center and grows with time t. Hunter (a) solution 
for t < describes an "exploding collapse" in which the mass of the central collapsing region gets smaller and smaller, 
converging to zero as t reaches 0, being surrounded by an expanding envelope. The central density still blows up to 
a "singularity" at t = 0. This solution can also be extended beyond t = to a late-time solution in which the whole 
cloud expands away; the singularity disappears and the central density keeps decreasing. 

Maeda and Harada |2^| found a unique relevant (unstable) mode for the Hunter (a) solution but no relevant mode 
for the Larson-Penston solution. It was concluded that the former is a critical solution while the latter is an attractor 
solution. In the former case, the eigenvalue k of the unique relevant mode was evaluated to k ~ 9.4637 by solving the 
eigenvalue problem. Other members of Hunter's family have more than one relevant modes, which implies that these 
solutions are not critical solutions. The homogeneous free-fall solution suffers from kink instability 27]. The value of 
the critical exponent 7 is ~ 0.10567. Since we have seen that the parameter for the family of initial data sets is given 
by the temperature T, the scaling laws in the present system are given by 

oc \T-T*p (2.26) 

for the supercritical case, and 

Pmaxcx \ T -T*\- 2 \ (2-27) 

for the subcritical case. It should be noted that the parameter is not necessarily the temperature but may be the 
mass or the radius of the cloud if we fix the temperature. 



III. NUMERICAL SIMULATIONS 



We have performed numerical simulations of the spherical collapse of an isothermal Newtonian gas. The basic 
equations for the evolution of the system are given by Eqs. (|2.1|l - (|2.4|l . The code is based on the Lagrangian description 
of hydrodynamics and a finite difference scheme. We have chosen Lagrangian (comoving) radial coordinates, attached 
to spherical shells with infinitesimal thickness. The mass inside coordinate f is consequently constant throughout the 
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simulation. We call such coordinates mass coordinates. Shells are initially positioned at constant intershell distance 
and their mass is computed from the desired density profile (we do not use equal mass shells). With these choices, we 
achieve high accuracy for the dynamics in the dense central region where the shells collapse during the simulations. 
The number of shells used in the simulation is 10 4 . One simulation with 10 5 shells has been run to check the influence 
of numerical resolution. No effect has been detected. This is proof of the accuracy of the code. In this section, wc 
adopt units in which G = k/m = 1 and the total mass mtot and initial radius r sur f,i of the cloud are unity. 

We have computed the evolution of four types of initial density profiles. Details of the models for the initial density 
profiles are summarized in Table [I] The initial velocity profile is always set to be identically zero. This assumption 
seems the most relevant for astrophysical applications, but the same dynamics should develop for different initial 
velocity profiles. Model 2 is our reference model, in which the initial density profile is given by a quartic polynomial 
and both the density and the density gradient (i.e., pressure forces) are continuous at the cloud surface. We have 
computed collapse simulations for various values of the temperature. (Actually, we have computed for various sound 
speeds and calculated the temperature because the sound speed is rather dynamically meaningful.) 

For each model, we find a critical temperature T* , which is of the same order of magnitude as the virial temperature 
estimated from the initial gravitational energy of the system. For temperatures just below T* , the cloud collapses 
and the central density goes to infinity in a finite time, producing a singularity. We call this collapse a supercritical 
collapse. For temperatures just above the critical temperature, the cloud begins to collapse but bounces afterward, 
and the density finally goes to zero everywhere. We call this incomplete "collapse" a subcritical collapse. The results 
for the critical temperture for each model are summarized in Table [I] In this table, we quote only the number of 
digits that are common between the critical values computed using the two different resolutions, 10 4 and 3 x 10 4 . 

In the critical and supercritical cases, we can estimate the time t s - mg of singularity formation at the center in 
different ways. First we can just define it as the instant when the numerical density reaches a very high threshold 
(e.g., p ~ 10 15 while initially of the order of one). Or we can assume the convergence to a self-similar behavior and 
use the relation p c oc (i s ing — t)~ 2 for the central density p c . Both approaches have been tried and give the same 
results. When the temperature is fine-tuned to the critical one, a case that we call the critical case for convenience, 
the simulations show that the collapse approaches the Hunter (a) solution for a while before it deviates. The finer 
the tuning of the temperature, the closer the approach and the later the deviation. Figure ^ shows density profiles 
in terms of the dimensionless quantity 47r / o(i S i ng — t) 2 for a tuning of the temperature to the critical temperature 
\6T\/T* = \T — T*\/T* ~ I0~ 15 , along with the theoretical Hunter (a) solution. In this case we estimate i s i ng using 
the p c oc (i s i n g — t) relation. This estimation is local in time and it gives the value which t s i ng would take if Hunter 
(a) solution was fully realized. We can see that the density profile converges to Hunter (a) solution. A discrepancy 
appears at large radii because the numerical configuration has a finite mass while the theoretical Hunter (a) solution 
has an infinite mass. This discrepancy moves to larger and larger value of r/(~c s t) as we approach t s - mg . Therefore, 
we associate the Hunter (a) solution to the critical case. 

In the near-critical supercritical case, the collapse first approaches the critical solution. Figure shows the evolution 
of the nondimensional central density parameter 

Q = ln[47rGp c (i sing -t) 2 ], (3.1) 

as a function of the central density p c itself. Here again, we estimate i S i ng locally. We compute i S mg,n at each time 
step n, using the following relation 

Pen (^sing.n ^n) — Pc,n— 1 (^sing.n t-n— l) ? C*^"^) 

where t n and p c n are the time and central density at the n-th step. We can then deduce Q. In this figure, the theoretical 
values for self-similar solutions, such as the homogeneous free-fall, Larson-Penston and Hunter (a) solutions, are also 
plotted. We find that the curve of the critical collapse approaches the value of the Hunter (a) (critical) solution. 
However, it afterward deviates and approaches the value of the Larson-Penston solution. This happens because the 
tuning of the temperature to the critical temperature is not perfect, and because the critical solution is unstable 
against such fluctuations. During and after the approach to the critical solution, the surrounding region expands 
away and leaves a finite-mass "core" in the central region. The density field in the central collapsed core at late times 
is well described by the Larson-Penston solution. This is the case early on for a supercritical collapse also, where the 
temperature is not fine-tuned to the critical temperature. This type of collapse also produce a finite (more massive) 
core and an expanding shell. We can see the convergence to Larson-Penston solution in Fig. [3J Therefore we identify 
the Larson-Penston solution as an attractor solution. 

Figure 0] shows the velocity profile around the center for a supercritical collapse. We can see that a finite-size 
collapsing core is formed. Moreover the physical radius of the collapsing core is almost constant. Figure [S] shows 
the time evolution of the collapsing mass, for different (supercritical) values of the temperature. We can see that at 
all supercritical temperature, the collapsing mass tends to a constant and therefore the "core" mass is well defined 
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at the time of the singularity. Figure [H] shows the relation between the core mass and the temperature fluctuation 
\T— T*\, The numerical results reproduce the scaling law derived by the renormalization group analysis in the previous 
section. The numerical value of the critical exponent is ~ 0.11, which agrees well with the theoretical value ~ 0.10567 
determined from the eigenvalue problem of the linear mode analysis. It should be noted that the power-law relation 
holds approximately for \ST\/T* < 0.01. 

For the subcritical collapse, the central density reaches a maximum and then bounces away. Figure shows the 
relation between the maximum density and the temperature excess T —T*. It shows that the maximum density also 
follows a scaling law, as is expected from the renormalization group analysis. The measured value ~ —0.22 of the 
critical exponent agrees with the theoretical value ~ —0.21134. Again, it should be noted that the power-law relation 
holds approximately for ST/T* < 0.01. 

Although we have concentrated on model 2, the above features have been checked for models 1 and 3. They are 
identical in all three models except for the specific value of the critical temperature. Model 4, the top-hat model, is 
special because the cloud is initially homogeneous. In this case, the central region of the cloud remains homogeneous 
and freely falls to a singularity if the condition TOtot/?"surf,i > 7r 2 Cg/2 is satisfied j2£j- If the condition is not satisfied, 
the collapse converges to Larson-Penston (attractor) solution. 



IV. ASTROPHYSICAL IMPLICATION 



In this section, we study the mass spectrum of collapsed cores resulting from temperature fluctuations around the 
critical case, according to the treatment of Yokoyama [llj. We build a toy model which shows how the critical gravi- 
tational collapse may affect the properties of the interstellar medium, and in particular the processes of fragmentation 
in molecular clouds and star formation. It should be emphasized that the purpose of this section is not to derive the 
whole shape of the observed initial mass function (IMF) but to show how a link can be established between an initial 
fluctuation field, and the resulting mass spectrum for the collapsed objects. Numerous additional physical processes 
would have to be introduced to derive an IMF. 

We will use 

M corc (T) = KM tot (l - ^ , (4.1) 

for the masses of cores formed by collapsing isothermal gas clouds with mass M tot and a temperature which is slightly 
less than T*. K is a constant and 7 ss 0.11 is independent of initial density profile. This relation is theoretically valid 
for I T - T* |< T*. We can check on Fig. ©that it is actually valid for \T - T*|/T* < 0.01 or M core < 0.3M tot . 
We assume a Gaussian probability distribution for the temperature fluctuations for simplicity, 



where T m is the temperature with the maximum probability and a is the dispersion, which has the dimension of 
a temperature. For simplicity, we consider a collection of clouds of the same mass and size with a distribution of 
temperature obeying the above probability function. Equation (|4.2|) allows us to compute the fraction of clouds 
collapsing to cores 

(3= [ P(T)dT. (4.3) 

JT'il-K- 1 /-!) 

The lower integration limit reflects the fact that a core with a mass heavier than the initial gas cloud cannot be 
formed. 

We can now determine the mass sprectrum of cores formed from a collection of gas clouds with masses M to t- We 
define the mass spectrum as the number dN of cores per logarithmic mass bin, normalized by 



InMtot J^/- 



This mass function is given by 
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dN 



d(lnM core ) 



-r'pfTfMco,,)) 

T* ( M corc 




r 1 - 



M 



KM t , 



i/ 7 ' 



(4.5) 



From numerical simulations, we have obtained 7 ~ 0.11 and K ~ 1.2. There are two parameters p = T*/T m 
and g = cr/T,„ in the mass spectrum. The mass spectra are depicted for (p, q) — (1.0,0.1) and (p, q) — (1.0,0.7) in 
Fig. Eta), while for (p,q) = (0.1,0.1) and (p, q) = (0.1,0.7) in Fig.EIb). Core masses heavier than the total mass are 
cut off. The mass M p associated with the peak in the spectra is smaller than M COTe = M to t for (p, q) = (1.0, 0.1) and 
(0.1, 0.1), while equal to M COYe = M to t for (p, q) = (1.0, 0.7) and (0.1, 0.7) [these regions are however not fully relevant 
because of the restrictions to be applied to Eq. I|4.1j . M p and the average mass M av of the formed core defined as 



C t0t *** 



dN 

■ rfA/cor 



■dM r , 



r 

Jo 



M t , 



dN 
dM caI , 
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(4.6) 



are summarized in Table [H] Although M p can be just equal to M to t, the average mass M av is always smaller. 

The main feature of the mass spectrum is a steep slope in the low-mass tail, due to the small value of the critical 
exponent 7. This fact has several implications. First, the cores are created with nearly the mass of the initial gas 
cloud. In the formalism originally introduced by Press and Schechter |30). it is assumed that overdense regions with 
a single mass scale, if they collapse, result in bound objects with the same mass. However, our result implies that the 
mass of the collapsed object is not the same as the mass of the initial overdense region, but may be rather smaller if the 
collapse is near-critical. Second, the critical collapse alone cannot account for the shape of the observed IMF, which 
shows a negative exponent on most of the mass spectum while we find a large positive exponent. Other phenomenon 
are responsible for the negative exponent. 



V. SUMMARY 



We have analyzed the critical behavior in the collapse of a Newtonian isothermal gas by the renormalization group 
approach and in numerical simulations. The renormalization group analysis in the Newtonian case is similar to the 
general relativity case. We have observed the critical phenomenon in the numerical simulation, which can be compared 
to those observed during black hole formation in the collapse of radiation fluid in general relativity. The scaling law 
for the formed black hole mass with critical exponent ~ 0.36 in general relativity is replaced by a scaling law for 
the collapsed core mass with critical exponent ~ 0.11. The critical solution is a continuous self-similar solution, 
specifically the first member of Hunter's family of solutions. 

For the supercritical collapse, the collapse near the center converges (early for low temperature, late for almost crit- 
ical temperature) to another continuous self-similar solution, the Larson-Penston solution. Generic collapse converges 
to this attractor solution, which suggests the universality of self-similarity in gravitational collapse. This phenomenon 
is an example for the validity of the self-similarity hypothesis in gravitational collapse. 

Some implications of Newtonian critical phenomena to astrophysics have been discussed. The critical behavior 
alone does not determine the IMF, clumps mass spectrum in molecular clouds and so on. However it suggest some 
limitations to Press-Schechter's model. The isothermal gas is a good approximation for a gas system with a short 
cooling time, while it will not be valid for optically thick clouds, where adiabatic gas model are more relevant. 
Consequently, critical phenomena should be investigated for more general equations of state. Research along this line 
is now under investigation |31| . 

Newtonian gravity is a simple model for more complicated gravitational theories, such as general relativity. New- 
tonian critical behaviors in gravitational collapse provide a good laboratory for investigating various effects on the 
critical phenomena in general relativity, such as deviation from spherical symmetry, rotation, external environment, 
a variety of matter fields, and so on. 
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TABLE I: Models of initial density profiles. 



Model 


Density Profile^ pi 


Critical Temperature T* 


f 


oc (f - r 2 ) 


0.4882608 


2 


cc (1 -r 2 ) 2 


0.5351886 


3 


cx cos 2 (7rr/2) 


0.5470739 


4 


<x6(l-r) 


0.2752568 



t The total mass and initial size of the cloud are normalized to be unity. 



TABLE II: M p and M av for the formed cores. 





M av /M tot 


Mp/Mtot 


(1.0,0.1) 


0.87 


0.94 


(1.0,0.7) 


0.90 


1.0 


(0.1,0.1) 


0.86 


0.95 


(0.1,0.7) 


0.90 


1.0 



10 












FIG. 1: Snapshots of density profile are plotted in terms of the dimensionless quantity AnGpt 2 as a function of r/(—c B t), for the 
critical value of the temperature in model 2 (\5T\/T* = \T — T"* | /r" ~ 10~ lj ) at different times. The origin of time coordinate 
is set to be the time of singularity formation. The theoretical Hunter (a) solution is also plotted. Numerical profiles converge 
to the Hunter (a) solution. 
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Hunter A 




Central density 

FIG. 2: Evolution of Q = ln[47rGp c (ismg — t) ] as a function of the central density p c , for model 2. The time tsing of singularity 
formation is estimated locally. The theoretical values of Q for several self-similar solutions are also plotted in this figure. See 
text for details. 
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FIG. 3: Snapshots of density profile are plotted in terms of the dimensionless quantity 4nGpt 2 as a function of r/(— Cgt), for 
a supercritical value of the temperature in model 2 (8T/T* — —1/2) at different times. The origin of time coordinate is set to 
be the time of singularity formation. The theoretical Larson-Penston solution is also plotted. Numerical profiles converge to 
the Larson-Penston solution. 




FIG. 4: The evolution of the velocity profile for a supercritical collapse (model 2), (a) in linear scale and (b) in logarithmic 
scale. The temperature is about 1% less than the critical value. A collapsed core forms in the center, surrounded by an 
expanding shell. Label t\ is the time when the central density reaches 10, and label ti is the time when we have stopped the 
calculation. 
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(b) 

FIG. 5: Evolution of the collapsing mass as a function of time for different values of supercritical temperatures [(ST = T — T*) 
difference to the critical temperature] for model 2, (a) in linear scale and (b) in logarithmic scale. The collapsing mass is defined 
as the total mass of shells having a negative radial velocity. For all values of the temperature the mass reaches a constant value 
early on, forming a well defined "core" . 
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FIG. 6: Scaling law for the mass of the collapsed core in supercritical collapses. The masses of the collapsed core formed in 
supercritical collapses in each model are plotted for different temperatures just below the critical temperature. The theoretical 
slope determined by the renormalization group analysis is also depicted. 
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FIG. 7: Scaling law for the maximum density reached in subcritical collapses. The maximum densities reached in subcritical 
collapses in each model are plotted for different temperatures just above the critical one. The theoretical slope determined by 
the renormalization group analysis is also depicted. 
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FIG. 8: Mass functions with four sets of parameter values for p = T*/T m and q = cr/T m for 7 = 0.11 and K — 1.2. In (a), the 
solid and the dashed lines correspond to (p,q) = (1.0,0.1) and (p,q) = (1.0,0.7), respectively. In (b), the solid and the dashed 
lines correspond to (p, q) = (0.1,0.1) and (p, q) = (0.1,0.7), respectively. 



